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ABSTRACT 


The  two-dimensional  indentation  of  an  elastic  half  space  by  a 
rigid  punch  under  a  slowly-applied  normal  load  is  considered,  for 
the  case  in  which  there  is  a  finite  coefficient  of  friction  p  between 
the  surfaces.  The  contact  area  is  then  divided  into  an  inner  ad¬ 
hesive  region  -c  <  x  <  c  in  which  the  surface  displacements  are 
known,  surrounded  by  regions  c  <  |x|  <1  in  which  the  friction  is 
limiting  and  the  lateral  displacement  (which  must  increase  in 
proportion  to  the  overall  load)  is  not  known  in  advance.  The 
problem  is  formulated  in  terms  of  a  coupled  pair  of  singular  integral 
equations  for  the  normal  and  shear  stresses  <r  and  t  at  the  surface; 
these  are  combined  to  give  a  single  homogeneous  Fredholm  equation 
with  positive  kernel  for  a  quantity  4*  proportional  to  the  difference 
t  -  \ur  in  the  adhesive  region.  The  largest  eigenvalue  of  this 
equation,  for  which  <$>  >  0,  gives  the  adhesive  boundary  c  in  terms 
of  p  and  Poisson's  ration  v  .  A  similarity  transformation  shows  that 
c  has  the  same  value  for  both  flat -faced  and  power  law  punches. 


AN  EIGENVALUE  PROBLEM  FOR  ELASTIC  CONTACT 
WITH  FINITE  FRICTION 

D.  A.  Spence^ 


1.  Introduction 

The  problem  posed  by  the  Indentation  of  an  elastic  half  space  by  a 
rigid  punch  has  a  straightforward  solution  in  the  case  of  frictionless  contact, 
for  which  a  boundary  condition  of  zero  shear  stress  beneath  the  punch  applies. 
As  a  step  towards  the  physical  situation  in  which  the  shear  stress  is  con¬ 
trolled  by  friction,  solutions  have  also  been  found  in  a  number  of  cases  for 
the  opposite  limit  of  fully  adhesive  contact,  in  which  no  relative  slip  between 
the  mdentor  and  the  half- space  is  allowed  once  contact  has  been  established. 
The  simplest  such  case  is  the  symmetrical  two-dimensional  indentation  of  a  half 
space  by  a  flat  punch  over  the  interval  -l<x  <1  considered  by  Muskhelishvili  (1953 
pp.  475-7).  The  solution  of  the  linear  elastic  equations  in  this  case,  although 
physically  realistic  near  the  centre  of  the  punch,  is  not  so  near  the  edges, 
where  the  ratio  of  shear  stress  t  to  normal  stress  <r  at  the  surface,  given  by 

where  k  is  a  material  constant,  is  divergent  and  oscil¬ 
latory.  This  behaviour  would  be  modified  either  by  slip  or,  if  this  were  pre¬ 
vented,  by  plastic  deformation  near  the  edges. 

In  the  present  paper,  the  same  problem  is  treated  under  the  more  general 

assumption  of  finite  friction,  with  coefficient  p,  between  the  surfaces.  The 
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case  of  the  flat  punch  is  indicated  in  figure  la  (p.  33).  We  suppose  that  over  a 
central  part  C'C  of  the  contact  area,  the  friction  is  sufficient  to  prevent  slip 
taking  place  as  the  load  P  is  slowly  increased,  but  that  in  the  remaining 
parts  A'C'  and  CA,  the  friction  is  limiting  so  that  the  material  is  displaced 
inwards  as  well  as  downwards  under  progressive  application  of  the  load.  It 
is  assumed  that 

(i)  the  load  is  applied  sufficiently  slowly  for  static  equilibrium  to  apply 
at  each  stage; 

(li)  it  is  applied  monotonically,  so  that  slip  is  always  in  the  same 
direction. 

At  the  outset,  it  is  not  clear  that  a  solution  of  the  simple  form  indicated 
will  exist,  since  there  might  have  been  more  thar  one  region  of  adhesion, 
with  interspersed  regions  of  slip  (and  this  may  well  be  the  configuration  to  be 
looked  for  in  the  corresponding  "unloading"  problem).  But  we  are  able  to 
show  that  the  problem  as  posed  possesses  a  unique  solution  satisfying  the 
linear  elastic  equations  and  the  physical  boundary  conditions.  The  first  of 
these,  over  the  interval  C'C,  is  imposed  in  the  form  of  an  inequality  on  the 
stresses 

M  <  jxlo- 1  (1.  1) 

together  with  a  requirement  of  zero  displacement  parallel  to  the  interface, 
while  on  the  intervals  A'C'  and  CA  we  set 

I T I  =  Jilo-I  (1.2) 

but  cannot  specify  the  parallel  displacement  u,  which  must  finally  be  found 


as  part  of  the  solution. 
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p  is  then  effectively  the  limiting  value  of  the  coefficient  of  sliding 
friction  as  the  relative  speed  goes  to  zero.  This  may  possibly  be  lower  than 
the  coefficient  of  static  friction,  but  if  the  inequality  (1.1)  is  satisfied  when 
p  has  its  sliding  value,  it  necessarily  holds  also  if  the  static  coefficient 
is  higher. 

The  extent  of  the  slip  region,  given  by  the  ratio  CC'/AA'  =  c  say  is 
an  eigenvalue  of  the  problem  which  cannot  be  specified  in  advance,  but  we 
find  that  for  given  values  of  the  physical  constants  (p  and  the  Poisson  ratio 
v)  there  is  just  one  value  of  c  for  which  the  inequality  can  be  satisfied, 
c  is  identified  as  an  eigenvalue  of  a  certain  homogeneous  Fredholm  equation 
with  a  positive  kernel,  and  the  general  theory  of  positive  linear  forms  shows 
that  only  one  eigenfanction  of  such  an  equation  is  everywhere  non -negative, 
so  c  can  be  uniquely  identified  as  the  corresponding  eigenvalue. 

Our  main  object  in  the  present  paper  is  to  formulate  the  mathematical 
problem  leading  to  this  eigenvalue,  and  to  determine  its  value  precisely.  This 
provides  an  analytical  check  on  the  accuracy  of  numerical  methods  which  are 
potentially  able  to  give  detailed  information  about  the  distributions  of  stress 
under  the  punch  more  easily.  For  this  reason  the  stress  distributions  are  not 
explicitly  calculated  here  except  in  limiting  cases.  A  numerical  approach  to 
the  problem  that  will  be  reported  separately  is  indicated  in  section  5.  3. 

The  governing  equations  are  set  out  in  section  2.  Before  proceeding 
to  the  solution  for  the  flat  punch,  we  note  also  the  form  taken  by  the  equa¬ 
tions  for  the  more  general  case  of  a  power-law  indentor  z  =  Brn  as  indicated 
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in  figure  lb  (e.  g.  a  circular  cylinder  (n  =  2)  or  a  flat  wedge  (n  =  1),  and 
show  by  appeal  to  a  similarity  argument  given  in  an  earlier  paper  (Spence 
1968a)  that  the  solution  for  the  stresses  in  this  case  can  be  obtained  by 
quadrature  from  that  for  the  flat  punch.  It  follows  that  the  extent  of  the  slip 
region  is  the  same  in  both  cases,  despite  the  different  forms  of  the  distribu¬ 
tion  of  normal  stress.  Thus  the  present  solution  of  the  eigenvalue  problem 
covers  the  whole  range  of  self-similar  indentors  with  symmetrical  loading. 

A  useful  limiting  approximation  to  the  slip  problem  can  be  obtained  by 
neglecting  the  change  in  the  normal  stress  from  its  frictionless  distribution 
due  to  the  presence  of  shear  forces.  This  is  calculated  in  section  2.  2  and 
leads  to  an  expression  for  c  similar  in  form  to  that  given  by  Galin  (1945)  in 
an  approximate  but  much  earlier  treatment  of  the  present  problem.  In  section  3, 
the  full  linear  equations  are  treated;  these  are  reduced  to  a  single  Fredholm 
equation  for  a  quantity  r(x)  proportional  to  the  difference  plo- 1  -  (t)  within  the 
adhesive  region,  treating  its  extent  as  known.  Achange  of  variable  relates  c  to  an 
eigenvalue  of  the  Fredholm  equation,  as  shown  in  section  4,  and  theresultsof  cal¬ 
culations  of  the  numerical  value  of  c  as  a  function  of  p  and  v  are  outlined  in 
section  5. 

2.  Governing  equations 

Suppose  ax,  ay  are  physical  coordinates  parallel  and  normal  to  the 
surface  of  an  elastic  half-space,  the  surface  values  of  the  normal  and  shear 
stresses  being  denoted  by 

(Vy=0=-(P/a)P(x)’  {Txy)y=0  =  (P/a)q(x) 
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where  P  is  the  total  applied  normal  force  over  the  contact  area  |x|  <  1,  so 
that 

1 

J  p(x)dx  =  1  (2. 1) 

-1 

and  let  au(x),  av(x)  be  the  corresponding  surface  displacements  due  to  the 
load.  Then  it  follows  from  results  of  Muskhelishvili  (1953,  chapter  19)  that 
p,  q,  u,  v  are  connected  by  the  coupled  pair  of  singular  integral  equations 


1  r  P(t)dt  .  ,  ,  ,  g  .  a  dv 

7  /  j  ToT  vq(x)  =  ■(f-V)p  dx 


(2.2) 


1  r  q(t)dt  ,  ,  ,  G  ,  a  du 

-/  ffr  -vpw=  (rrr> — 


p  dx 


(2.  3) 


in  which  y  =  (1  -  2v)/(2-2v),  v  being  Poisson's  ratio  and  G  the  shear 
modulus.  If  the  load  is  applied  symmetrically  about  x  =  0,  p  and  u'  are 

even  functions  of  x,  q  and  v'  are  odd. 

% 

The\?oundary  conditions  on  the  normal  displacements  are 


-1<x<1:  ( 


_S-v£  dv  = 
1-v  '  p  dx 


n-1 


(flat  punch,  figure  la)  (2.4) 

(sgnx)  (power  law  indentor  z  =  Br°  , 
figure  lb)  } 


where  A=  (r^-)( — )  nB  is, a  constant,  initially  unknown,  relating  the  geomef.y 
1-v  p  x 

V 

of  the  indentor  and  the  size  ofcthe  contact  region  to  the  normal  force. 

The  lateral  displacement  within  the  adhesive  region  is  zero  for  the  flat 
punch;  and  for  the  power  law  indentor  similarity  considerations  of  the  type 
given  for  the  fully  adhesive  indentor  by  Spence  (1968a)  continue  to  apply  in 
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the  case  of  partial  slip, 


t 


giving  the  second  boundary  condition 


-c  <x<c: 


( 


G 

1-v 


a  du 
p  dx 


(flat  punch) 
(power  law)  # 


(2.  5) 


In  the  remaining  intervals  c  <  |x|  <  1,  can  not  be  specified.  The 

dx 

values  of  A  and  C  are  eventually  to  be  determined  from  the  conditions  that 
the  stresses  vanish  at  the  boundaries  x  =  ±1  of  the  contact  region.  (For  the 
flat  punch,  by  contrast,  the  stresses  become  infinite  at  these  edges). 

The  ooundary  conditions  on  the  stresses  are  neatly  expressible  in  terms 
of  the  difference 


r(x)  =  p(x)  -  ^^(sgn  x)  .  (2.  6) 

f* 

To  prevent  slip  in  the  adhesive  region  we  require 

r(x)  >  0  (-c<x<c)  (2.7a) 

while  in  the  regions  of  slip 

r(x)  =  0  (c  <  |x |  <  1)  .  (2.7b) 

In  these  boundary  conditions,  c  is  an  eigenvalue  to  be  determined  in  terms 
of  the  physical  constants  jx,  y  . 


2.1.  Reduction  of  "power  law"  equations  to  the  flat  punch  problem 

Before  investigating  the  eigenvalue  problem,  we  note  that  the  equations 
for  the  power  law  indentor  are  reducible  to  those  for  the  flat  punch.  Using 
the  symmetry  of  p  and  antisymmetry  of  q  to  write  them  on  the  interval  (0, 1), 
t  This  point  will  be  discussed  further  in  a  subsequent  report. 
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they  become  for  x  >  0 


_  .  2x  r  P(t)dt  .  n-1 

0  <  x  <  1:  —  f  1  -  +  yq(x)  =  -Ax 

-1  t  -x 


(2.8) 


0  <  x  <  c:  ~  f  -  y  p(x)  =  Cx”'1 

0  t  -x 

and  the  condition  that  the  stresses  vanish  at  the  edge  is 

p(l)  =  0  =  q(l)  . 


<2.  9) 


(2.10) 


Application  of  the  differential  operator  “  to  these  equations  removes 

n  n  ax 

the  constants  A  and  C,  » educing  them  to 


1  p  (t)dt 

0<x<U  —  f  ■  °  2  *  [y<I0M]  =  # 

0  t  -x 


(2.11) 


2  1  *<J  Ct)dt 

0<X<C!  ^  /  -2— 2--VP0W  =  ° 
0  t  -  x 


respectively,  where 


(212) 


P0(x)  =  P(x)  -  ”(xp)'f  qQ(x)  =  q  -  ^(xq)1  . 


(2.13) 


These  are  exactly  the  equations  that  apply  in  the  flat  punch  case,  so  the 
solution  for  the  power  law  indentor  is  given  by  quadrature,  using  the  boundary 
condition  (2. 10),  as 


.  1  pn(t)dt  .  1  q  (t)dt 

,  ,  n-i  r  o  ,  .  n-i  c  o 

P(X)  =  nx  J  — - - ,  q(x)  =  nx  J  — - - 


(2.14) 


x  t 


x  t 


If  the  distribution  of  p^(x)  is  normalised  in  accordance  with  (2. 1),  that  of 
p(x)  is  automatically  so,  since  integration  of  (2.13)  gives 
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1  1 
/  p(x)dx  =  f  p.(x)dx  . 

0  0 

Moreover,  the  distribution  of  p  and  q  so  found  satisfy  the  frictional 
boundary  conditions  (2.7a  and  b)  with  the  same  value  of  c  as  for  the  flat 
punch,  since 


.  .  n-1 

r(x)  =  nx 


r  <  x  <  1 
0<x<c  . 


Therefore  the  solution  of  the  eigenvalue  problem  for  the  flat  punch  also 
gives  the  solution  for  a  power  law  lndentor,  and  we  shall  from  now  on  only 
treat  the  former  case,  and  work  from  equations  (2. 11),  and  (2. 12). 


2. 2.  A  limiting  solution 

Since  Iq/pl  is  at  most  equal  to  k,  the  term  in  square  brackets  in 
equation  (2.11)  is  of  order  py  times  the  normal  stress  term,  and  a  first  ap¬ 
proximation  that  throws  some  light  on  the  physics  can  be  found  by  omitting 
the  term,  when  (2. 11)  has  the  normalised  solution 


PQ(x)  -  (l/w)(l  -x2)'l  .  (2.15) 

This  approximation,  equivalent  to  neglecting  the  effect  of  shear  forces  on  the 
distribution  of  normal  pressure,  has  been  used  for  the  calculation  of  the  shear 
forces  themselves  In  other  contoxts  (e.  g.  Goodman  1962  ,  Mindlln  1949).  On 
substitution  of  (2.15),  (2. 12)  becomes  a  singular  integral  equotion  for  qp(x)  on 
the  interval  0  <  x  <  c,  since  we  can  write  qQ  ■  npfl  in  (c,  1): 


(2.16) 


i 


2  rCtVt)dt 

-  f  -j-r  *  *w  * 0 

0  t  -x 

where 

,  1  tp_  (t>dt 

-r— r-vpo(x>  • 

0  t  -  x 


Inversion  on  the  interval  (0,  c)  gives  the  solution  bounded  at  x  =  0  as 


v*>-T/l44)ijrv • 

0  x  c  - x  /  y-x 


(2.17)  I 


The  solution  is  in  general  unbounded  at  x  =  c,  but  is  bounded  if 
C  2  24 

/  (c  -  y  )  1  x  (y)dy  *  o,  and  since  one  of  the  terms  of  x  proportional  to 
0 

p  and  the  other  to  y,  the  vanishing  of  this  integral  gives  an  expression  for 

i  \ 

the  ratio  p/y: 


H/y  *  K(c)/K*(c), 


(2.18) 


*/2  2  2  -4  / - 2 

where  K(c)  is  the  complete  elliptic  integral  j*(l -c  sin  e)'*d0,  and  K'(c) »  Kh/l  .c  ) 

0 

Galln's  result,  obtained  by  a  lengthy  analysis  using  approximations  in 
a  complex  plane  that  are  difficult  to  interpret  physically,  I i  of  the  same  form. 

x 

In  his  equation  12.  36,  the  right  hand  side  is  Identical  with  the  above,  but  the 
left,  in  oxir  notation,  is  (tan"V)/i  which  is  the  same  when  p  and 

y  are  both  small  compared  with  1  .  The  accuracy  of  the  two  solutions  relati  /e 
to  the  exact  solution  is  indicated  in  figure  2.  Galln's  Is  very  close  for 

i 

v  «  0.  3,  i.  e.  y  =  2/7,  and  reasonably  so  for  v  =  0,  y  =  equation  (2. 18) 
is  not  as  close  as  Gatin's  but  is  still  useful  numerically,  and  knowledge  of 
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1 


the  limit  as  p-*  0,  y—0  with  \i/y  fixed  is  also  useful  in  the  subsequent 

1 

analysis. 

\ 

Asymptotic  expressions  when  c  is  small,  and  when  c  is  close  to  1  , 
are  found  from  the  limiting  behaviour  of  the  elliptic  integrals: 


K(c)/K'(c)  ~ 


(2/ir)  log(4/tyl-c  )  as  c  •*  1 


(w/2)  log(4/c) 


as  c  —  0 


so  that  the  width  of  the  adhesive  region  is  given,  according  to  (2. 18),  by 


4e_TT^/2^  as  p-*0  for  fixed  y 

1  -  as  p  -*  «  for  fixed  y  . 


(2.19) 


1 


When  (2.18)  is  satisfied,  a  calculation  using  elliptic  functions  (appendix 
D)  shows  that 

qQ(x)  =  [h/ttK(c)]  (1  -X2)"^  Ffsin'1  c)  (2.20) 


e 


2  2  - 1 

where  F(£,c)  is  the  incomplete  elliptic  integral  f  (1-c  sin  0)‘*de  .  In 

0 

the  further  limit  c  -*■  0,  F(|)  -*  £  and  K(c)  n/2,  so 


q0(X)  2  ,  -lx 

pp.(x)  "  17  8  n  C 


(2.21) 


in  the  adhesive  region  0  <  x  <  c  .  Figure  3  shows  computed  curves  of  q  /pp 
for  the  case  v  =  0  for  three  values  of  c,  for  comparison  with  this  expression. 
The  agreement  is  best  for  the  smallest  value  of  c  .  The  computed  curves  for  * 
0.  3  would  lie  approximately  mid  viray  between  those  for  y  =  0  and  the  broken 


curves. 
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3.  Reduction  of  the  dual  problem  to  a  Fredholm  equation 

In  this  section  we  treat  the  flat  punch  equations  (2. 11),  (2. 12)  but  will 
from  now  on  omit  the  suffix  0  .  If  q  Is  eliminated  In  favour  of  r,  these 
become 

0  <  x  <  It  p(x)  ♦  2x/*)(cot  w  o)f  «  r(x)  (3.1) 

0  y  -x 

0  <  x  <  ct  -y2p(x)  ♦  (2/*)(tanw  o)  /  ^ ydy  ■  0  ,  (j.2) 

0  y  -x 

where  tan  wa  =  (0  <  a  <  ))  , 


which  are  to  be  solved  subject  to  the  requirement  that  r  vanishes  on  (c,  1) 
and  Is  bounded  and  non -negative  on  (0,  c)  In  accordance  with  (2. 7).  A 
singularity  In  p  at  x  =  1  can  however  be  admitted  in  the  solution. 

In  the  physical  problem  It  would  be  natural  to  fix  the  values  of  a  and 
y,  and  try  to  find  c  .  However,  It  Is  more  convenient  analytically  to  pro¬ 
ceed  for  the  present  as  If  all  three  were  known,  but  we  shall  find  In  due 
course  that  a  solution  to  (3. 1),  (3.  2)  satisfying  the  further  requirement  that 
r(x)  is  bounded  and  non-nogatlve  on  (0,  c)  exists  only  when  r  certain  com¬ 
bination  of  a,  y  and  c  Is  equal  to  the  first  eigenvalue  of  the  Fredholm 
equation  (4.7)  below.  (3.1)  Is  a  singular  Integral  equation  for  p(x)  In  terms  of 


r(x),  with  solution  obtainable  by  standard  methods  (appendix  A)  as 


where 
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(An  arbitrary  multiple  of  l/xu(x)  forms  part  of  the  solution,  but  must  be  ex¬ 
cluded  since  it  is  unbounded  at  x  =  0  .  The  remaining  integral  is  regular  at 

2  -  —  4-(j 

x  =  0,  but  behaves  like  (1  -  x  )  2  near  x  =  1. ) 

This  expression  for  p(x)  can  now  be  substituted  into  (3.2)  to  yield  a  homog¬ 
eneous  equation  for  r  alone,  as  follows: 

(i)  On  substitution  of  p  into  the  second  term  of  (3.2),  the  equation 
becomes 


-y  p(x)  - 


(  sin2iro> 

|  r°  yr(y)dy  f 

l  -  J 

2  2  -4[ 
0  y  -  x 

v2  1  2 


(ii)  The  order  of  integration  on  the  right  hand  side  may  be  reversed  by 
use  of  the  Bertrand -Poincard  lemma  (Tricomi  1957  p.  172)  when  the  term  becomes 


(sin  it  a) 


-r(x)  +  /  w(u)  r(u)du  f 


2  , 

y  dy 


2  2  2  2 
<*>(y)(y  -x  )(u  -y  ) 


The  inner  integral  here  can  be  split  into  partial  fractions  and  evaluated  by 
use  of  a  general  result  derived  in  Appendix  B  : 

1 


2x 


2/  sin 


n  TTO  \  r 

*  )Jo 


_dy_ 


2  2 

w(y)(y  -x  ) 


xcostto’  2  ,  .  2, 

~Zw~-x  g(a)A(x  >’ 


(3.  5) 


where  A(x2)  =  F(l,  1+a;  x2)  and  g(o-)  =  r(  {+  a)/r(l  +  a)  r(|)  . 

(iii)  When  these  results  are  combined,  (3.4)  reduces  to 


.saas/r,  .*a}  w*  ♦  ig(. 

0  ^  s  u  -  x 


z  z 

)  siniro  J  G(x,u  )w(u)r(u)du, 
0 


where 


oo 


2  2%  r  2  24  2  2vl  2  2,  v  ,  2n  2n-2  2,  t  2n  , 

G(x  ,  u  )  =  [x  A(x  )  -  u  A(u  )]/ (x  -  u  )  =  2^  a  (x  +x  u  +  .  . .  +  u  ),  (3.  6) 

n=0  n 
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and 


lc 

a^  =  (  $),.  /(I  ♦  e)^  Is  the  coefficient  of  z  In  the  expansion  of  A(z)  . 
The  first  term  and  the  first  part  of  the  second  term  combine  to  produce  exactly 
-p(x),  by  (3.  3),  while  the  remainder  of  the  second  term  cancels  the  second 
term  on  the  left  of  (3.4),  which  therefore  reduces  to 

2  2  ^  2  2 

0  <  x  <  c:  (1  •  y  )p(x)  -  -  g(a)(slnwa)  f  G(x  ,  u  )  w(u)  r(u)du  =  0  (3.7) 

*  0 

Since  p(x)  has  already  been  expressed  in  terms  of  r(x),  (3.7)  is  equivalent 
to  a  Fredholm  equation  for  r  on  the  Interval  (0,  c)  .  It  is  reduced  to 
canonical  form  in  the  next  section. 

In  die  slip  region  c  <  x  <  1,  the  right  hand  side  of  (3.7)  equals 
(  ^^-)(  p^5x  ’  displacement  u(x)  could  therefore  be  calculated 

once  the  equation  has  been  solved  for  r(x)  .  It  can  be  confirmed  directly  from 
(2.  3)  that  ^  and  therefore  u(x)  are  negative,  i.  e.  that  the  displacements  in 
the  slipping  regions  are  inwards.  This  is  consistent  with  the  positive  sign 
of  q  for  x  >  0  (shear  stress  away  from  the  origin)  since  the  frictional  force 
opposes  the  motion. 

4.  The  Fredholm  equation  for  r(x) 

In  (3.  3)  and  (3.7)  write 

X2  «  c2(,  (1  -x2)^'®  r(x)  *  ♦(£),  (l  -  x)^*  °  p(x)  •  *(»)  ,  (4.1) 

Then  the  equations  become 

+({) .  «<(>  .^o....i"..)y,ili-.ibjda  M  ,2) 

■  14  »«y » 
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and 

Mi)  «  gi°i,lnwo  cZe  f  V*'1  H(i,i)>4(nXli,  (4.J) 

*(1  -  Y  i  0 

in  which  c  plays  the  role  of  an  eigenvalue,  to  be  fixed  so  that  (4. 2)  and 
(4.  3)  possess  a  solution  such  that 

0  <♦<*)<•  (4.4) 


on  (0, 1).  The  kernel  In  (4.  3)  Is 


H(e,i))«a-c2e)*‘*G(c2ttc2n)«  £  c2(l>*n,h  e“nn  <«•*> 

m,  n«0 


say. 


The  inverse  of  (4.  2)  that  Is  bounded  on  (0, 1)  is 


(4.6) 


[This  excludes  an  arbitrary  multiple  of  fc^/0-t)^"**  to  comply  with  (4. 4)J. 
(4.  3)  and  (4.  6)  give 


*Mi)  ■  /  n*‘l  *U,n)4(n>dn, 


(4  7) 


where  K<t,  r\)  • 


and 


(^)L-',{,H,{.„.(^f)  Ecl'''V«On"  ,4 


mn 


H  .  (1  -  Yl(/c**(s(.)lJ  . 


14  91 


Here 


4w(t)  •  L'l(tm)  -  Im40<l>. 


(4  10) 


.  m -l 

♦  g(#)(cot**>tMl*C)  %  ,r  |lf 


r»0 
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where  4>Q(£)  =  L  V)>  which  is  evaluated  using  the  result  in  Appendix  B,  and 
can  be  written  in  two  ways: 


♦0<t>  =  [r(l-a)/r(i)r(|-a)]^(l  -i)^'“F(l,l-a;|-a;l-e) 
=  1  -  2[r(l-a)/r(i)r(i-a-)]£^  F(i,  . 


(4.11) 


The  first  shows  that  <j>  and  therefore  all  the  d>  tend  to  zero  like  (l-£)2 
as  £-*  1;  the  second  shows  that  =  1,  while  4>m(0)  =  0  for  all  m  >  0  . 

Insertion  of  these  expressions  in  (4.  8)  shows  that  K  is  bounded  on  0<£,r|<l  , 
and  if  we  write  $(£)  =  <\>  (£)  the  equation  is 


M>*(£)  =  f1  (§Ti)"2(1"0f)K(|,  T1)<t>*(T1)dTi  (4.12) 

0 

in  which,  because  of  the  boundedness  of  K,  the  kernel  is  square  integrable 
provided  a  >  0,  so  the  standard  Fredholm  theory  shows  that  solutions  exist 
for  at  most  a  discrete  set  of  eigenvalues  X  >  \  >  \  .  .  .  say.  Moreover, 

1  ”  u  J 

all  the  coefficients  h  are  >0  (this  is  confirmed  by  a  detailed  calculation 

mn 

in  §5),  and  ^^(^J  >  0  on  the  open  interval,  so  K  is  positive,  and  is 
the  Perron  root,  which  is  positive  and  is  the  maximum  in  absolute  value;  and 
the  corresponding  eigenfunction  $^(£)  (i.  e.  the  solution  of  (4.7)  when 
X  =  Xj)  is  the  only  eigenfunction  which  is  positive  on  the  whole  interval 
(0,  1)  .  This  therefore  gives  the  solution  of  the  physical  problem.  Since  the 
kernel  K  depends  on  c  and  a,  we  find  X^  (c,  a)  from  (4.7)  and  must  finally 
determine  c(a',y)  numerically  by  inserting  this  value  for  X  in  (4.  9).  Asa 
first  step  we  replace  the  integral  equation  by  a  set  of  linear  algebraic  equations. 
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Solution  of  the  Integral  equation  (4.7) 


Hie  separable  form  of  the  kernel  suggests  a  solution  of  the  form 


♦*»■£  «m  c-ve.  • 


m=0 


(4.13) 


Then  since  the  { +  }  are  linearly  Independent,  the  coefficients  of  cme  (£) 

m  m 

on  the  two  sides  of  the  equation  can  be  equated.  This  gives  a  set  of  linear 

equations  for  the  {a  }  which  can  be  written 

m 


\a  « y 

m 


n,  P=0 


m+2n.  .  p 

:  h  1  c  o 

mn  np  p 


where 


sinira  t ^  e+n-1 


1  s 

np  *g(a)  J  '  p 


0 


(4.14) 


Now  a  detailed  calculation  In  appendix  C  shows  that 


1  =  h  , 

np  np  ’ 


(4.15) 


so  (4.14)  can  *  <  written  in  matrix  form  as 


*8 


(4.16) 


where  K  =  CHC,  H  *  {h  },  C  *  [  c  ,  |  being  infinite  matrices. 

mn  I  _c 


Since  K  and  K2  are  positive,  If  is  the  Perron  root  of  K2  It  follows 
that  »JT^  Is  the  Perron  root  of  K,  the  corresponding  eigenvector  a ^  and 


eigenfunction  ♦  1  (£)  *  Y e  1  ♦  (£)  being  the  same  in  each  case.  Moreover, 

m  m 

+^(()  Is  the  first  eigenfunction  of  the  clmpler  equation 


.16. 
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(4.17) 


=  /  V*’1  K(t,i))+(vixifi 
0 


where  -c^)*1  «  £  c2m  *  (e)r,m  . 

'  '  m=0 


Because  of  (4. 15),  it  follows  from  (4.  8)  that 


/  K<e,U  «  K(6,t,)  . 

0 

[An  equivalent  result  Is  J  -c2fc  t)*1  K(&,  rj)dt  =  H(4>  n)> 

0 


(4.18) 


(4.19)] 


So  all  eigenfunctions  of  (4.17)  are  also  eigenfunctions  of  (4.7).  The  reverse 
does  not  follow  In  general  although  it  is  true  for  the  Perron  root.  It  may  be 
noted  that  (4. 17)  can  also  be  written  in  the  form 


i  a-i 


(4.20) 


in  which  the  left  hand  side  equals  *7  X  «(,(£),  but  the  author  has  not  been  able 
to  derive  this  equation  directly  from  the  physical  problem. 

5.  Calculation  of  the  Perron  root  Kj  and  the  physical  eigenvalue 


The  matrix  H  formed  by  the  coefficient  h  (m,  n  =  0,  1,  2, 

mn 

equation  (4.  5)  is  the  product 


. )  in 


H  =  BA 


(5.1) 


1  a^  a2  ... 


where  A  =  I  al  a2  a3 


a2a3  a4  •• 


bl  1 
b2  bl  1 
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are  respectively  an  infinite  Hankel  matrix  and  an  Infinite  lower  triangle,  the 
a's  being  the  coefficients  In  the  expansion  of  A(x)  already  quoted  (equation 
(3.  6J)  and  the  b's  those  In  the  expansion  of  (1  -xf  'a,  namely 

bk  *  •  (5-2) 


The  fact  that  all  elements  of  H  are  positive  Is  deduced  as  follows: 

(1/  The  m**1  element  In  the  left  hand  column  is 

h  *  b  +  b  .a.  *  ...  a  =  coefficient  of  xm  In  expansion  of  (1  -x)^’aA(x)  . 
mo  m  m  - 1  i  m 

This  function  equals  F(|  +  a,  a;  1  ♦  a:  x)  by  Kummer's  transformation  (Erdllyl 
etal  1953,  p.  105  $2.9  (2)), 


whence 


_<l_  lLlsIzl  >  o 

m m!  ’ 


(5.  3) 


(11)  h  ,  .  -  h  =  -b  a  >  0  for  all  m  >  0  (since  all  the  b’s  ex- 

m-i,  n+i  mn  m  n 

cept  hQ  are  negative)  so  the  elements  of  H  Increase  to  the  right  along  lines 
of  constant  m  +  n,  so  are  all  positive. 

The  Perron  root  of  K  Is  ^  kj,  and  may  be  computed  by  the  power  method 
of  successive  multiplication  by  K  until  convergence  in  achieved.  In  numerical 
work,  K  was  first  truncated  to  containing  the  first  N  rows  and  columns, 
leading  to  a  value  the  procedure  was  then  repeated  for  N  ♦  5  until 

the  value  of  ^/k^  did  not  change  by  more  than  10 .  With  increasing  c  , 
the  value  of  N  to  secure  convergence  Increased  from  10  at  c2  ■  0.1  to  45 
at  c2  -  0. 95. 

Finally,  the  solution  of  the  physical  problem  in  the  form  of  e  as  a 
function  of  c  and  y  Is  found  by  Iteration  from  (4.9),  written  logarithmically 


In  the  form 
.18- 
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or  = 


log  c  f 


.-JLlogtf  --Y 

log  ^(c,  a)  ~ 


(5.4) 


2a 


+  -  logg(a) 
a 


Starting  from  a  trial  value  of  a  on  the  right,  with  c  and  y  held  fixed,  this 
converges  very  rapidly.  Results  are  given  in  table  1. 

The  friction  coefficient  is  then  found  as 

P  s  (1/y)  tan  Wflf  • 


Values  of  p  so  found  from  the  values  of  a  given  in  table  1  are  listed  in 
table  2  for  the  cases  *  =  0  and  v  a  o.  3,  and  plotted  as  solid  lines  in 
figure  2  . 


5. 1.  Limiting  form  of  solution  as  a  —  0 

log  Vj  * 

In  the  limit  a  -*  0,  ■  -  ■  tends  to  a  limit  x  (c)  say,  as  can  be  in  - 

ferred  directly  from  the  matrix  (see  footnote  page  21),  or  from  the  limiting 
solution  of  section  2.2. 

The  limit  can  be  approached  by  allowing  y  to  tend  to  zero,  and  since 
(1/a)  logg  •*  -  log  4,  (5. 4)  gives 


lim  1  -  \  j (logc  +  p* -log4)  . 

Y-0  [_log(l  -  y  )J  / 

The  limit  on  the  left  is  p/y  =  K(c)/K'(c)  by  (2. 18),  whence 


P*(c)  =  logi  -  j  K'(o)/K(c)  . 


(5.5) 


(5.6) 
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A  convenient  table  is  given  on  page  608  of  Abramowitz  and  Segun. 

Table  1  "hows  that  (log  Xj)/6r  differs  appreciably  from  \  only  when 
c  is  close  to  1  and  y  close  to  £;  for  v  =  0.  3  it  is  sufficiently  accurate 
to  write  *  in  place  of  (log  \j)/ar  so  the  determination  of  from  the  matrix 
can  be  avoided,  although  Iteration  to  find  a  is  still  necessary  since 
(logg)A  s  .log  4  is  not  sufficiently  accurate.  The  full  expansion  is 

00 

logg(ar)  *  -a log 4  +  £  (-2 a)"  ^  (5.7) 

n=2 

where  t](n)  is  the  sum  of  reciprocal  rfi1  powers  with  alternating  signs 

2  JL 

(rj(2)  =  ir  /12),  and  if  one  extra  term  is  retained,  (5. 1)  can  then  be  written 

.^ogd-vV^pgi-f]'1  (5.8) 

which  leads  to  values  of  p  in  4 -figure  agreement  with  those  tabulated  in 
most  cases.  This  expression  was  used  to  calculate  the  entries  in  table  2 

i 

for  c  <  0. 7  . 

The  solution  for  $(£)  can  be  written  down  explicitly  when  a  -*•  0,  since 
in  this  limit  the  first  eigenvector  of  the  matrix  K  is  (1,  0,  0,  ... )  (This  is 
because  in  the  limit  all  the  elements  of  the  left  hand  column  of  K  are  zero, 
except  for  k0Q). 

Therefore,  as  a  —  0 

♦(1)(e)  =  inn  ye) 

or-*  0 

=  1  -  |  F(i,  1;  3/2;  i\  by  (4.11)  , 

II 

=  1  -  -  sin"1  (5.9) 

ir  * 
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in  perfect  agreement  with  (2.  20)  when  allowance  is  made  for  the  change 
of  notation. 


Footnote; 

To  infer  the  existence  of  the  limit  directly  from  the  matrix,  note  that  all 
the  elements  of  the  left  hand  column  of  K  except  the  first  contain  a  as  a 
factor.  Therefore  if  we  write 

=  1  +  2a  X  , 

as  «-*0,  X  tends  to  X  (c),  given  by 


where  the  stars  denote  the  limits  of  the  quantities  as  a  -*  0,  so  ( h  /<* )a  = 

(A)  /m.  m!  etc. 
m 

I  have  not  found  a  direct  derivation  of  (5.  6)  from  this  determinant,  but 

$ 

have  satisfied  myself  that  the  two  expressions  for  X  are  the  same  with  6 
figure  accuracy  by  means  of  a  computation,  which  took  only  a  fraction  of  a 
second  of  machine  time. 
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5.  2.  Asymptotic  expression  as  c  —  1 


The  writer  has  also  carried  out  an  analysis  parallel  to  that  reported  here, 
in  which  the  problem  is  cast  into  the  form  of  an  integral  equation  for  the  unknown 
strain  u'(x)  in  the  slip  region  e  <x<l  .  This  is  tractable  by  the  Wiener  - 
Hopf  technique  when  c  is  close  to  1,  and  as  in  the  present  approach,  an 
eigenvalue  problem  for  c  follows  from  the  requirement  of  bounded  stress.es 
at  x  =  c,  giving  rise  to  the  condition 

-i  . 

=  tanh  o  +  *  (a,  k)  +  log  2  (5.10) 

2 

with  an  error  of  order  (1-c)  ,  where 

k=  (2/n)  tanh  =  (1  /n)  ln(3  -4V) 

and  *  (a,  k)  =  (1/k)  Im  log[r(l  -  a  -  -y  )/r(l  )] 

2  2,2 
=  -jy  [1  +  0.  2326  irar+  —  (a--— )  +  ...]  . 

Values  of  p  calculated  from  this  expression  are  included  in  table  2  and  are 
seen  to  be  extremely  accurate  for  c  greater  than  about  0.6. 

5.3.  Numerical  solution  of  the  coupled  equations 

Mention  may  also  be  made  of  a  numerical  attack  on  the  problem,  which 

will  be  described  separately  in  an  MRC  Technical  Summary  Report  by  the 

x  x 

author.  If  we  write  y  for  x,  and  apply  the  operators  J  >— — —  ,  J  i— — — 

0  \Jx  -y  0  v  x  -y1* 

to  (2.11)  and  (2.12)  respectively,  they  are  put  in  the  form  of  coupled  Volterra 
equations : 

1  tpQ(t)dt 

0  <  x  <  1:  I  ■  1  ■  +  y 

x  <7  tZ  -  x2 


/ 


xyqQ(y)dy 
0  -y2 


(5.11) 
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1 


(5.12) 


1 


which  are  particularly  suitable  for  numerical'  solution,  as  the  earlier  singular 
equations  are  not.  They  have  been  solved  numerically  on  the  Unlvac  1108 
computer  at  the  University  of  Wisconsin  Computing  Center,  using  piecewise 
constant  interpolations  to  p,  q  with  up  to  64  sub -intervals,  the  condition 
q  *  |ip  being  app’ied  on  the  portion  c  <  x  <  1  with  an  initially  guessed  value 

i 

of  p  .  For  each  fixed  c,  an  iteration  was  performed  until  the  value  of  I1 

I 

was  such  that  q  -  pp  was  positive  throughout  (0,  c)  and  tended  to  zero  at 

\ 

c  .  In  this  way  the  polrits  shown  by  circles  in  figure  2,  and  the  corresponding 
stress  distributions  in  figures  3  and  4  were  obtained. 

t  am  very  grateful  to  Professor  Ben  Noble  for  suggesting  the  limit  looked 
at  in  section  2.2  and  for  his  continued  interest,  land  to  Verlyn  Erickson  for 
programming  the  calculations. 

i 


! 
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Carlemann  (1922)  showed  that  tha  general  solution  of  the  equation 


a(s)  u(s)  -  X  f  -M?*-  e  g(t) 
-1  * 


(X  >  0)  where  a(s)  and  g(s)  are  analytic  for  -1  <  *  <  1  can  be  written 

u(>)  r  a(g)9(»)  ,  Ke T(S)  f 1  e'T(t)g(t)dt  KeT(8)  ; 

b2(s)  b(S)  -1  b(t)  (t  •  8)  1  -  8  ’ 

where  K  Is  an  arbitrary  constant,  and 


b2(s)  =  a2(s)  +  (Xir)2  . 


Equations  (2. 16),  (3. 1)  and  (4.  2)  are  all  of  this  form,  and  the  inverses  quoted 
are  obtained  by  suitable  changes  of  variable  and  choice  of  K:  eg.  in  (3.1) 
write  x2  =  £(s  +1),  y2  =  }(t+l),  x_1p(x)  =  u(s),  x-1r(x)  =  -g(s),  when  the 
equation  becomes  identical  with  (Al)  with  a  =  -1,  wX  =  cot  va  .  Then  e(s)  = 
|+o,  expr(s)  =  +  b(s)  =  cosec  ira  and  the  solution  (A2),  translated 

back  into  the  original  variables,  is 


(cosec  iro)p(x) 


2°  yr(fldy 

2  2 
y  -x 


-2a.,  2  -i+a 

+  jKx  (1  -  x  )  2 


If  we  choose 
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the  last  two  term*  combine  to  give  the  expression  for  p<x)  quoted  in  (J  S|# 


which  Is  regular  at  xiO, 


Appendix  B.  Evaluation  of  the  integrals  (3.  5)  and  (4. 11) 

These  are  special  cases  of  the  following  general  result  for  real  a,  b  <1, 
0  <  s  <  1: 


costra 
s  (1-s) 


simra 


/ 


1 


dt 


m-b) 


0  ta(l-t)b(t-s) 


F(l,a+b;a+l;s) .  (Bl) 
(0-) 


r(l+a)r(l-a-b) 

1  'V_/  Q  U 

This  is  obtained  by  evaluating  the  loop  integral  -r—r  f  dz/z  (z-1)  (z-s) 

Zvl 

+  oo 

round  the  two  sides  of  the  positive  real  axis  in  the  z  plane,  with  semi-circular 
indents  above  and  below  the  point  z  =  s  .  The  second  term  on  the  left  hand 
side  is  the  contribution  from  the  two  sides  of  the  interval  (0, 1),  and  the  first 
that  from  the  indents  at  z  =  s  .  The  right  hand  side  ccmes  from  the  integral  on 
(1,  oo),  which  can  be  transformed  into  Euler's  integral  for  the  hypergeometric 
function.  (In  the  special  case  a  +  b  =  0,  the  contour  can  be  replaced  by  a 
closed  loop  surrounding  the  points  0,  1,  and  the  right  hand  side,  which  re¬ 
duces  to  1,  is  the  residue  at  infinity. )  From  this  result 

2  2 

(i)  (3.  5)  is  obtained  immediately  by  writing  a  =  a,  b  =  j  -  a,  s  =  x,  t  =  y, 

(Li)  To  derive  (4.11),  write  41  =  1  in  (4.  6)  whence 

1 


♦  (e)a2(H.)2 

V9'  sin  ™ 


sin  na 

"1 - r 

e2a-i); 


cos  tt a 
t r 


/ 


“P" 


0  t1*(1-t1)*'£V£) 


Use  of  (Bl)  with  a  =  \  -  or,  b  =  j,  s  =  1  -  £,  t  =  1  -  r)  shows  that  the  quantity  in 
square  brackets  is 


[r(i)/r(a)r(|  -  a)]  F(l,  1-a,  |  -a;l-£)  , 


from  which  the  first  line  of  (4.11)  follows  since  r(a)r(l-a)  =  ir/simra  . 

The  second  line  is  obtained  from  the  first  by  writing  a  =  0,  b  =  a,  c=£ 
in  lines  (22),  (43),  (1)  and  (17)  of  Erd^lyi  §2.9,  PP.  105-108. 
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Appendix  C.  Evaluation  of  1  _  siriTr^  r  aip-1 

pq  Trg(a)  J  *  YqVb  b 


Substitution  from  (4.10)  and  the  first  line  of  (4.11)  gives 

1 

(* 

0 


oo 


i 


pq 


/  |p+q+a- *(!-!)*■"  £ 


(l-a)n  (1.|)nd€ 


T(i+o)r(j-a) 


1  i  i  „  q-i 

f  tH  i-e)1  Z 


r=0 


n=0  (f  -a) 

2  n 


r+cH-p-1 

a  t  dt  . 

q-r-1  * 


(1) 


(2) 


On  termwise  integration, 


,  oo  (1-a) 

(1)  =  ar(p  +  q  +  or+^  £  n 


r(£  +  a)  r(p+q+2)  n=0  (P+<*+2)n 


The  sum  of  the  series  is  F(l,  1  -a;  p+q+2;l),  which  equals  (p+q+l)/(p+q+a) 
by  Erdelyi  1953,  p.  61(14),  so 


i  p+q 

m  =  <—£ _ )  (a+  ll£ta  -  Y  a  b 

D+d  +  O-  fD+dll  L>  as  b 


'p+q+o-'  (p+q)! 


s=0 


s  p+q - s 


by  (5.3) 


Similarly,  integration  of  the  second  term  gives 

[-1 

C_ 

q 


(2)  = 


COSTTQ  -TT 

tt  q-r-1  cosn 

r=0 


q-1 


)  b  =  -)  a  b 
ira  p+r+1  s  p+q-s  . 

s=0 


Adding  the  two  expressions  gives 


pq 


p+q 

=  Z  a 

s=q 


„  b 

s  p+q. 


=  h 

s  pq 


by  (5.1)  . 


An  alternative  demonstration  of  this  result  is  possible  by  direct  integration  of 
various  hypergeometric  functions  to  obtain  the  equivalent  equation  (4.19),  but 
the  analysis  is  more  lengthy. 
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lx  D  -  Dllptlc  Integrals  for  limiting  stress  distribution 


If  the  order  of  Integration  It  reverted  in  the  flrtt  term  on  the  right  of  (2. 17  ^ 


we  heve 


q0(«i 


c  ^ 

•  7TT7\' Jjrrv*-''-'  P?-rV  • 

*  vc  -x  I  c  v  1  -  t  vl-y  y  -x  J 

A 


,  C  ,  2  2  1  .  .-1.2  2.! 

where  I  (x,t)  •  :|-L  when  0  <  x  <  c  <  t  <  1 

9  o  (t  -y2)(y  -x2)  r  -  x2 


V° ' 


2x  j  /■*.  t2-c2.J  dt  fC  c2-y2  1  dy 

L  c  l-t  t  -X  0  l-y*  /•«  J 


exhibiting  e  singularity  at  x  ■  c  unless  the  sum  Inside  the  square  bracket 
vanishes  at  this  point,  l.e.  unless 

r*  dt  r°  dy 

0  ■  ’»)  - 1 - r  ♦  v  J  - r - r 

c  (l-t2)  (t-c,  0  (l-y  )  (c  -y2) 

which  gives  p/y  ■  K(c)/K'(c),  equation(2. 18).  When  this  is  satisfied,  the 
previous  expression  can  be  rearranged  as 


q0(x)  c  pfj  JcZ-xZ  [pll(x>  -  y  I2 <x) ]  , 


where  Ij(x)  =  f  - r— ^ — 


r1  dt  1  ..  ,1-c2  I,  _2\ 

I  T  T  =  2  2*  **  1_C 

c  (l-t2,  (t2-c2)  (t2-x2)  1-x  *-x 


L  2  2  . 27 
(c  -X  )(1  -X  ) 


[E’F($,c)  +  K'E(e,c)  -  K'F(£,c)] 
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C  -  1  c2 

a  Jn  2w  2  2»#  2  2.  2  ni{  2  ’  C) 

0  V  cl -y  )(X  -y  )(y  -X  )  x  x 


and  IjW  ■  /  I—] 2 - ^ 


x*/(c2-x2)(l  -X2) 


(KE<4,c)  -  Er<e,c)|  . 


Here  If j  is  tho  comploto  elliptic  intogial  of  the  third  kind  (Byrd  end  Friedmann, 
*  1  x 

p.  225),  4  «  sin  -,  K  *  K(c)  and  E  =  E(c)  arc  the  complcto  integrals  of 

C 

the  first  and  second  kinds,  and  dashes  denote  the  same  integrals  with 

argument  Vl-c2  .  E(£,c),  F(4,c)  are  the  incomplete  elliptic  Integrals 

f  (1-c2  sln20)*i  de  .  Then  since  y  e  nK'/K,  these  results  combine  to  give 
0 

q„(x)  = - ^ - [KE'  ♦  K‘E  -  KK'l 

0  2  /— T  K(c) 

w  Vl-X 

The  quantity  In  the  square  brackets  equals  n/2  by  Legendre's  formula,  so 
qQ(x)  reduces  to  the  expression  quoted  (2. 19). 
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TABLE  1 


Iteration  to  find  a, 

for  v  >  0 

(V  o  J)  and 

V  a  0.  3(v  *  -  ) 

V  • 

0 

V  a 

0.  3 

c* 

N 

X*(C) 

'logx, 

1 

'logXj 

(equation  5. 5) 

2a 

a 

TT“ 

a 

.1 

10 

. 026080 

. 025997 

.059399 

. 026069 

.017  140 

.2 

10 

. 054620 

.054370 

. 070821 

.054586 

.020233 

.  3 

10 

.086190 

.085680 

. 080634 

.086155 

. 022824 

.4 

15 

.  121604 

. 120699 

. 090332 

.  121521 

.025303 

.  5 

20 

.  162072 

. 160555 

. 100748 

.  161932 

. 027886 

.6 

25 

.209533 

. 207058 

. 112806 

.209308 

.030765 

.7 

30 

.  267427 

.263323 

. 128027 

.  267064 

. 034229 

.75 

30 

. 302242 

.296837 

.  137751 

.  301768 

.036342 

.8 

35 

. 342929 

. 335596 

. 149923 

. 342307 

. 038884 

.85 

35 

. 392360 

.381957 

. 166332 

.  391447 

.042121 

.9 

40 

.456536 

.440404 

. 191387 

.455181 

. 046695 

.95 

45 

. 552638 

. 521275 

.241681 

.  549429 

.054643 
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TABLE  2 


Friction  coefficient  u  as  e  function  of  c 


v  m  o  v  ■  o.  3 


1  s 

1*  1 

1  H 

c 

(equation  5.4) 

(equation  5. 10) 

(equation  5.4) 

(equation  5. 10) 

.0$ 

. 20952 

.  10728 

.  1 

.25064 

.  12771 

.  15 

.28356 

.  14384 

.2 

. 31313 

.  15815 

.25 

.  34118 

.  17154 

.  3 

. 36869 

.  18449 

.35 

. 39633 

.  19733 

.4 

.42466 

.21028 

.45 

.45422 

.22356 

.5 

.48557 

.23738 

.55 

. 51941 

•  25199 

.6 

. 55660 

.26766 

.6325 

. 58332 

.27882 

.7071 

. 65504 

.  65722 

. 3074 1 

. 30819 

.7746 

.74003 

.74142 

. 33934 

.  33973 

.8367 

. 85080 

.85168 

. 37782 

.  37805 

.8660 

. 92393 

.92465 

.40135 

.40150 

.8944 

1.01845 

1.01914 

. 42969 

.42978 

.9220 

1.  15190 

1.  15260 

.46587 

.46597 

.9487 

1.  37198 

1.  37300 

.  51715 

.  51720 

.9747 

1.  89811 

1.  90518 

.  60680 

.  607  56 
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^  |  x*c  A(x=l) 

Slip  Adhesion 


Adhesion 


Figure  1:  Frictional  indentation  by  rigid  punch 
(Schematic) 


1 


8 


1.0 
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mal  stress  distributions  for 
-dimensional  flat  punch 


